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It is a well-known fact that, among the electroweak corrections, QED radiation 
gives the largest contribution and the needed precision requires a re-summation of the 
large logarithms which show up in perturbation theory. For annihilation processes, 
e + e~ — > f f, initial state radiation is a definable, gauge- invariant, concept and one has 
general tools to deal with it; the structure function approach and also the parton-shower 
method. However, when one tries to apply the algorithm to four-fermion processes that 
include non-annihilation channels a problem is faced: is it still possible to include QED 
corrections by making use of the standard tools? A systematization is attempted of 
several, recently proposed, algorithms. In particular, it is shown that starting from the 
exponentiation of soft photons one can still derive a description of QED radiation in 
terms of structure functions, i.e. the kernel for the hard scattering is convoluted with 
generalized structure functions where each of them is no longer function of one scale. 
Each external, charged, fermion leg brings a factor where a is the fine-structure 

constant, < x < 1 and A is a function which depends on the momenta of the charged 
particles. 



PACS Classification: 11.10.-z; 12.20.-m; 11.15.Bt; 12.38.Bx; 02.90. +p 
*Work supported by the European Union under contract HPRN-CT-2000-00149. 



1 Introduction 

In QED, and more generally in the standard model (SM), the most important terms in radiative 
corrections to various processes at high energy or at large momentum transfer are those that 
contain large logarithms of the type In (Q 2 /m 2 ), where Q 2 stands for a large kinematical variable. 
The quantity m is the mass of a light, charged particle, which emits photons, for example, the 
leptons or the light quarks. The origin of these large logarithms is connected to the presence of 
collinear divergences in the perturbative expansion. For QED in particular they arise whenever 
a photon is radiated in the direction of an incoming (outgoing) light fermion and we compute 
an exclusive observable. Collinear divergences are present both in virtual and in real corrections. 
As far as virtual corrections are concerned, they appear as a consequence of a mass singular 
renormalization procedure. Suppose we have a Green function F (q, . . .), where q is some external 
momentum of a loop diagram. Then the renormalization procedure is mass singular if 

F ren (g,...) = F(g 2 ,...)-F(0,...), (1) 

that is, the subtraction has been performed at q 2 = 0. 

The collinear divergences enter F Ten through F (0, . . .), in the limit — q 2 ^> m 2 , where 

-q 2 

F ren ~ln n -V , n=l,2. (2) 

When we consider, instead, the real corrections, terms like A -1 = (p — k) 2 + m 2 will arise as a 
consequence of a fermion of momentum p and mass m emitting a photon of momentum k. The 
singularity will arise if we set m = 0. With p 2 = k 2 = and k \\ p the denominator in A goes to 
zero. If we now integrate over k, then the expression for the diagram in which A appears becomes 
singular. Strictly speaking, collinear divergences are only present in theories where massless quanta 
couple to each other, which is not the case for QED, or for the SM, where the theory becomes 
almost massless, at high energies; that is, almost mass singular, due to the presence of low mass 
fermions. In those cases where we need a very high accuracy of theoretical predictions, the presence 
of large logarithms calls for a re-summation of the perturbation series, in spite of the low value of 
the coupling constant. The method for this summation was actually developed within QCD and 
is based on the factorization theorems, which allow us to split the contributions of large and small 
distances. As a result of this re-summation procedure one starts with the cross-section for the 
so-called hard process; that is, the process with large kinematic variables, which is subsequently 
convoluted with the structure function of the initial (final) particles JIJ]. 

For a given scale Q, the hard part of the cross-section is determined by distances of the order 
1/Q and is expressed in terms of the coupling constant at these distances, a{Q 2 ). This hard 
cross-section contains no large logarithms and is evaluated in perturbation theory. The whole 
contribution of large distances is inserted into the structure functions, which obey an equation of 
the renormalization-group type. 
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The structure functions describe the radiative corrections to a given process in the so-called 
leading logarithmic approximation (LLA): 



a a Q 2 

- L = - In — , 3 

7T 7T 7Tl Z 

where Q represents some large scale inherent to the problem under consideration. Therefore, in 
LLA we re-sum all terms of the form (a/irL) n . Sometimes it is possible to go beyond the LLA to 
guarantee the correct evaluation of the next-to-leading terms, a/ir (a/nL) n . The latter, usually 
included in the so-called .ff -factor, are not universal but require instead a comparison with the 
explicit calculation of the cross-section up to the two-loop level. 

However, for a general process we are not always in this favorable situation. Very often, 
therefore, one can find the statement that the choice of the appropriate scale in the structure 
functions is mandatory. This is a jargon for 'implementing the correct exponentiation factor in 
multi-photon emission'. 

The outline of the paper will be as follows: in Section |] we introduce the problem of a general 
treatment of QED radiation for arbitrary processes which do not allow a gauge invariant defini- 
tion of initial state radiation. In Section |3.1| - [3.3| we review the classic Yennie, Frautschi and 
Suura exponentiation algorithm (hereafter YFS), completely reviewed in terms of dimensional 
regularization for both virtual and real infrared terms. In particular we give the derivation of 
the YFS virtual form-factor in Section [|. In Section [5] we shown that the YFS real form-factor 
can be exactly evaluated in n-dimensions without having to split a soft part from the hard one. 
In Section |] we introduce an approximation to the exact result which is much simpler to handle 
in practical computations. The extension from two emitters to an arbitrary number is discussed 
in Section [7]. As a consequence of our approximation the result will be naturally expressed in 
terms of generalized structure functions with an overall exponent which is discussed in Section |j. 
Additional refinements are introduced and discussed in Section |]. Some general considerations 



are introduced in Section |10| that introduce the problems that one encounters in going beyond the 
lowest order result. Conclusions are shown in Section [U]. Finally, the most relevant integrals used 
in the paper are explicitly shown in appendix. 



2 Background of the problem 

The great success of high precision LEP physics is intimately linked to the possibility of splitting 
initial state QED radiation from the rest in a meaningful way and of using accurate determina- 
tions of the structure functions for the incoming e + e~ beams, the so-called s-channel structure 
functions g. 

For LEP 2 physics, although the required theoretical precision is not comparable to what 
we need around the Z resonance, we are not in the same fortunate situation O. First of all, 
initial versus final state radiation is not longer a meaningful concept. Secondly, a large fraction 
of processes are not dominated by annihilation and, therefore, the standard methods of using 
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s-channel structure functions fail to reproduce the correct result. Nevertheless, the language 
of structure functions is a useful one and it is desirable to include at least the bulk of large 
radiative corrections. Hence, structure functions are still applied for these processes, but some 
large uncertainty remains, connected with what is usually referred to as the problem of selecting 
the right scale. In a word, the choice of the energy scale is not a trivial issue. 

For processes where some exact perturbative calculation exists the scale inherent to structure 
functions can be determined by matching the two languages but this is not possible in general. The 
keyword in all these cases is, according to commonly accepted jargon, to 'select a suitable scale 
without knowing the exact one (two) loop calculation'. Needless to say this attempt is utopistic, 
although some ingenuous strategy has been devised in recent times. In particular we refer to some 
interesting work that can be found in ref. || and in ref. || 

More or less, all these attempts amount to start with some sort of naive soft + virtual photon 
approximation, not the one commonly employed in the Yennie, Frautschi and Suura formalism ||, 
but rather something as in the following expression: 

where m, {jpj) are the mass(momentum) of z-th charged particle, AE is the maximum energy of 
the soft photon (the boundary between soft- and hard-photons), E is the beam energy, and Qi 
the electric charge in unit of the e + charge. The factor €j is —1 for the initial particles and +1 for 
the final particles. The indices i, j run over all the charged particles in the initial and final states. 

We will say that any formulas as in Eq.(|4]) includes all double-logarithms. This terminology 
is not ordinarily accepted, so we stress that it means including, prior to integration, all terms 
proportional to 

\n(AE/E) \n(Q 2 /m 2 f ), (5) 

where Q 2 denotes, generically, some large scale present in the problem. It must be stressed that 
Eq.(|j) misses to incorporate virtual corrections and that the single-logarithmic part is omitted. 
The adopted strategy is to implement a structure function approach on the Born kernel for 
the process with the opinion that one is able to make an educate guess about the scale in these 
structure functions, something as in 

o"tot(s) = J dx_dy^dx + dx u dx d D e -(x_,Q 2 _)D e -(y_,Q 2 F _) 

x D e +(x+, Q 2 F+ )D u (x u , Ql)D d (x d , Q 2 d )a (s), (6) 

where all scales, Q 2 _ etc are guessed. Note that, for definiteness, we have taken a specific example. 

Another problem arises in adapting structure functions to t-channel processes. A non-accelerated 
charged particle cannot radiate. Consider now an incoming electron of momentum p- m , an outgoing 



da soft (s) da (s) 



dQ 



dQ 



exp 
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electron of momentum p out and the corresponding radiation of any number of photons from both 
legs, before and after the hard scattering, with total momentum K. Furthermore, let Q = Pi n —p ou t 
and consider processes dominated by the region where Q 2 ~ 0. The adopted strategy would be to 
fold the hard scattering cross-section with structure functions at virtualness Q 2 , 



where x is the probability of finding the electron within an electron with longitudinal momentum 
fraction x; however, the limit Q 2 — > (no radiation) will not be correctly reproduced. Of course, 
one can switch off the radiation for Q 2 < Qmin with Q^m chosen ad hoc, but this represents an 
artificial solution. Another solution is to replace the logarithm of the structure function with a 
power law behavior for small values of Q 2 but a continuous solution valid in any regime and not 
only for Q 2 m 2 e or Q 2 <C ml is desirable. 

Furthermore, all approaches that are missing virtual photonic corrections - at least the uni- 
versal, process independent, YFS form-factor - simulate their effect by imposing an effective lower 
energy cutoff on the photon energy, E c , and require ha.(E max / E c ) = 0(1), where E max is the 
maximum photon energy. 

The situation is slightly better for those t-channel processes where a perturbative calculation 
exists, see for instance the two-photon process e + e~ —>■ e + e~fi + fi~. Here the comparison between 
structure functions and the O (a) calculation |7| allows us to select a scale t once a ^-factor is 
included, box diagrams are neglected (but luckily enough their contribution is strongly suppressed) 
and one stays away from the forward scattering region where large deviations are expected and 
seen. 

Our approach is aimed to systematize this basic idea and, therefore, it must be clearly stated 
that it does not represent an attempt to move beyond the correct treatment of those terms that 
may lead to double-logarithmic enhancement. Simply, we start from the fundamental process of 
exponentiation of Yennie, Frautschi and Suura and try to understand its correct interplay with 
the language of structure functions. Complete virtual corrections are not included, therefore only 
the universal YFS factor is included, collinear single-logarithms are again missing since, basically, 
only soft photons enter into the scheme. However, the issue of 'selecting the scale' turns out to be 
much less arbitrary than in any previous approach. 

In particular, there is a situation where the exclusion of hard photons represents a bad approx- 
imation to the exact result. Consider two final state emitters, then owing to the Kinoshita, Lee 
and Nauenberg || theorem (hereafter KLN) the corresponding corrections are always small and 
free of large logarithms. They can become more sizeable only for tight cuts on the invariant mass 
of the emitting pair. The typical example is in the corrections to e + e~ — > // + //~ where the O (a) 
inclusive corrections are represented by the well-known factor 3/4 (a/V). This result, however, 
requires adding the matrix element for single hard bremsstrahlung and the complete O (a) virtual 
corrections. For a general process with many fermions in the final state the knowledge of this 
matrix element is usually missing and, unless tight cuts are imposed on the invariant masses of 
the final state pairs, we cannot reproduce the correct KLN limit. 
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In principle hard photons can be included for arbitrary processes in the so-called collinear ap- 
proximation, namely hard photons are allowed only within a small cone surrounding each charged 
external fermion. In this case, however, we are not allowed to integrate over the whole phase 
space of the photon since, in this case, there is no gauge invariant leading - sub-leading splitting 
of collinear radiation and the latter, in any case, will not be suppressed. Even more, the complete 
O (a) virtual corrections are needed as well and this can only be derived on a process-by-process 
basis. 

YFS exponentiation has been widely used in the past || and we have no pretension to be 
adding any substantial improvement. The only goal of this paper is to clarify the extraction of 
structure functions from the YFS program, without having to introduce a soft-hard separation 
in the YFS form factor, and extending their validity beyond the asymptotic region where all the 
invariants of the process are much larger than all fermion masses. Therefore, the main difference 
is that we adopt a slightly different variant of the YFS-approach where no separation is made 
between the soft and the hard region; rather we re-formulate in modern language an old proposal 
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3 QED corrections 

The material in this section is well known |J and we go through it with the main motivation 
of establishing notations and conventions. As an example consider the process 

e + ( P i) + e-(p 2 )^J( Pi ) + f(p 3 ), (8) 

3.1 Soft photon exponentiation 

The cross section for the emission of an extra soft photon of momentum k is 

/ x 2 

a f d?k I ^ 6iPi 



4 tt 2 J k \p? A Pi-kj 
where dao is the non-radiative differential cross-section. The variables Q\ are defined as follows: 

in — part. 9 2 —Q e 

in — antipart. 9 1 +Q e / in \ 
out — part. #3 +Qf 

out — antipart. # 4 —Qf 

and they satisfy conservation of charge, i.e. J2i @i — 0. We define the usual eikonal factor, 

Pi - K 
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which satisfies current conservation, 

j ■ k = °i = 0- (12) 

i 

Next we consider the generalization of process Eq.(|8]), e + e~ annihilation into several fermion- 
antifermion pairs and an arbitrary number of photons: 

e + ( P+ ) e-(p_) -> Ji + /i + + 7j + // + "7- (13) 

In our approach intermediate vector boson are unstable particles and never appear in the final state. 
When the emitted photons are soft the corresponding amplitude is approximated, by standard 
methods, and reads as follows: 

n 

ilC ft = (ie)"M J] <ka)-j(ka), 

a=l 

5M =E«< "HP (1 4 ) 

i=l,2l P* ' Ka 

where e is the photon polarization vector, e(k) ■ k = 0, and M is the non-radiative amplitude for 
the process. Note that, in soft approximation, only photons emitted by external charged fermions 
are relevant. This fact is connected with gauge invariance as it can be illustrated by considering 
e + e~ — > ZZ'-f. There are two diagrams in Born approximation, direct and crossed conversions. If 
we consider only emission from the external e ± lines, the Ward identity k ■ M ext = is violated. 
However k ■ M|° t ft = 0, so that M|° t ft is gauge invariant and the full identity is restore by including 
M int , i-e. emission from the internal electron line which, however, is neither infrared nor collinear 
divergent. 

The cross-section for Eq. (fH|) , for an infinite number of emitted soft photons, is 



i= E 1 [ dPS no ^ I3d dPS Iad J2 \ M n\ 2 (27rrV|p + +p_-]>>-X>;,], (15) 



21 

E 

n=0,oo "" " spins \ i=l j=l 

21 



dPS non - iad = (2 7if l J] d% 5 + (fi + 

i=l 

n 

dPS vad = (2 7r) 3n Y[d%5 + (k*), 5 + (k 2 ) = 6(k )5(k 2 ), (16) 
i=i 

where we have split the total phase-space into radiation and non-radiation phase spaces. We easily 
derive 

oo i . n 

E - / dPS non ^ d E I M o | 2 (2vr) 3 " J] d%5 + (k 2 ) 

n=0 U - J spins j=l 

(21 n \ 

p++p--Eft-E*i ■ (1 7 ) 
i=i 3=1 1 



According to the classical treatment one writes 



54 ( K ~ E k ij = / ^4 eX v{ iK ' X-t^ky x}, 



2/ 



K=p ++P _-J2q l . (18) 



i=l 

This result can be cast into the following form: 

21 



o ~ J fl d% 5 + (g 2 + mf) (2 tt) 6 '" 4 £ | M | 2 £ (p + + p_ - £ ft) 

i=l spins V i J 



(19) 



spins 

where the spectral function for the photon has been introduced: 



/n 
H d% 5 + (kf) I e ^(%) | 2 exp (- i kj ■ x) . (20) 

The photons can be re-summed to all orders giving 

r r/ 4 r 00 1 -in 

W = / ^p (iX ■ x) £ - [F(x) , 

l Z 71 J n=0 n - 

= — ^ y d 4 fc exp (-ik-x) 5 + {k 2 ) I e j M (A;) | 2 . (21) 

The Dirac delta-function, expressing four-momentum conservation, is therefore replaced inside 
Eq.(pO|) by the photon spectral function. The latter is defined through a Fourier-transform, 

E(K) = — / d 4 x exp (i K ■ x) E(x), 

(2 7Tj 4 J 

E(x)=exp{^ J d*k exp (-ik-x) 5 + {k 2 ) | f(k) | 2 }. (22) 

E(K) is the spectral function describing radiation and a = e 2 /47r is the fine-structure constant. 
An important property, following from charge conservation is 

and, therefore, we derive 

E(x) = exp{- ]T 0^- / d A k exp (-ik-x) 5 + (k 2 ) | - | 2 }. (24) 



Note that there is no delta-function expressing four- momentum conservation inside Eq.(^), not for 
the full process nor for the soft limit. Therefore, we are not allowed to use relations connected to 
momentum conservation. Later in this paper we will adopt a coplanar approximation for the exact 
spectral function, the difference between the two to be treated perturbatively, which will introduce 
again conservation but, this time, expressed in terms of the process where one incorporates photons 
emitted along the directions of the charged fermions. 



3.2 Re-organization of the perturbative expansion 



The exponentiation of infrared divergences is rigorous only in the limit where all photon mo- 
menta go to zero. The procedure adopted here will be to define an approximation to the exact 
cross section where the matrix element for the process are replaced by their soft limit. However, 
we do not introduce a cutoff separating hard from soft in the exponentiation, therefore the ap- 
proximated result receives contributions from photon momenta which are large as well as soft. 
Perturbation theory is then reorganized by evaluating the difference between the exact and the 
approximated contributions. The main goal of the present investigation will be to establish a link 
between our approximation and the structure function language. Note that, by virtue of Eq.fll2T), 
our approximation satisfies U(l) gauge invariance. 

Let M n denote the complete amplitude for the production of a final state {q} with the emission 
of n photons, 



M n (p+,P-;{q};k 1 



k n ) 



Pi- 



ll \M n 



(25) 



spin 



The entire procedure amounts to construct a perturbative expansion which starts with an ap- 
proximation that embodies the desired features of re-summation. We introduce, as usual, the 
factor 



J(k)=\ee(k)-j(k) 
Then the complete amplitude squared can be written as 

Pn 



Po n j ( k i) + e n j{ki)Pi{ki) +...+p n (h, ...,k n ). 

i=l l+i 



1=1 

A solution for the infrared-finite residuals (3 is 

A) = Po, Pi(k) = pi(k) -po J(k), etc. 
giving a perturbative expansion for the cross-section that starts with 

1 



(26) 



(27) 



(2? 



+ 



(2^) 4 
1 



(2tt) 



d 4 x exp (i K ■ x) E(x) J dPS q f3 
d 4 k5 + (k 2 ) exp(-ik-x) Pi(k) + ...]. 



(29) 
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Here dPS q is the phase-space for the final state fermions. At this point we can ask about the 
validity of a result that includes only f3 . The singular behavior of the amplitude for any radiative 
process is better examined in terms of the so-called dipole formalism. For simplicity let us consider 
the case of only two emitters, the generalization to an arbitrary number being straightforward, 
both outgoing. An approximation to the true amplitude squared that incorporates the correct 
singular behavior in the massless fermion limit is represented by 



M | 2 ~ | M non - rad | 2 gij 
9ij 



ypi-pj 



l-z(l-y) 



1 - z 



V= ^ r, z= Pi ' Pj (30) 

Pi ■ Pj + Pi ■ k + pj ■ k Pi-Pj+ Pj ■ k 

It is easily seen that the infrared limit corresponds to k — > or y — > 0, z —>■ 1 while the collinear 
limit to y — *■ 0, independently of z. Therefore the exponentiation procedure correctly accounts for 
photons that are infrared or infrared & collinear but not for photons that are hard & collinear. In 
the coplanar approximation, which receives contributions from photon momenta which are large 
as well as soft, we include all photons, soft or collinear that originates from 

Pi- Pi (31) 



Pi ■ kpj ■ k 

that, however, does not reproduce the correct results of Eq. (|30|) . In particular (3\ is essential to 
correctly reproduce the KLN result. 

3.3 Inclusion of virtual corrections 

Let Mo be the Born amplitude for the process where n photons are radiated. Let M\ be the 
same amplitude with one loop corrections included. Therefore we obtain 

M = exp (aB) \M + (Mj -aBM ) + ..], 

B = ~ Bi i ( 32 ) 

Bij is the universal, i.e. process independent, YFS virtual factor. Virtual corrections are opera- 
tively included through the following procedure: 

1) a — y\ exp (aB) | 2 a, 

2) Pq -> E sp ins [l M | 2 +2 ReM* (M 1 - a BM )] = f3 o + Pou ( 33 ) 

3) /3i -> p u 
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where the M in f3 denotes the one-loop corrected amplitude with no real photons and j3\ gives 
the Born amplitude with one emitted photon. The cross-section becomes 

a = J dPS q | exp(« B) \ 2 E (p + + p_ - ]T q) (/3 00 + (3 i) 

+ -—r I dPS q d 4 k5 + (k 2 )d 4 K \exp(aB) | 2 E(K) -J— f d 4 x 
(2 7r) 4 J " (2tt) J 

x exp{? (p + + p_ - J2 V - fc - K ) ■ x ) /M*) + ... (34) 



4 Evaluation of 5^ in dimensional regularization 



Consider an arbitrary process containing charged incoming and outgoing fermions, each with 
charge Qi and momentum pj. The YFS ^-function in Eq . (|32"D describes infrared divergent virtual 
photons associated with the external charged lines, therefore independent of the internal details 
of the process. Let us consider the definition of B^ in dimensional regularization, i.e. for 4: 



Bij = fi 



4— n 



k 2 



d n k r {2e iPi -k) tl {2t p 3 + k) 
k 2 ^k 2 -2eiPi-k k 2 + 2e j p j -k* 



(35) 



where e, = ±1. The connection between q and 0, is as follows: each particle has a charge Qi = 
0, —1, 2/3, —1/3, then 9i = e 4 Q iy or e« = ±1 for incoming(outgoing) fermion i (outgoing(incoming) 
anti-fermion i). Using 



k 2 - 2 6iPi ■ k = (k - €i p^ 2 + m- • 
k 2 + 2 6j pj ■ k = (k + 6j pj) 2 + m 2 



(2) , 

(3) , 



(36) 



we derive the following decomposition: 

1 



Bi 



+ 



(2) 2 
2 

MM 



4 (Pi-etPi-k) 



+ 



(3) ; 



4 [p] + e 3 - Pj • k) + A; 2 



4 e^e, pi ■ Pj + 2 (e* - Sj pj) ■ k - 



(37) 



which allows us to derive B^ in terms of standard one- loop scalar functions ||12|| . We have three 
terms 



B 



i.i 



5> 

i=i 



i 



where the first one can be written as 



4 = ^ 2 



d 



{B (p 2 ; 0, nm) - 2 [2p 2 B (p 2 ; 0, rm) + P- B x (p 2 ; 0, rm)] } 



(3? 



(39) 
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The various terms are computed as follows fi3|l : 



B o (-m 2 ;0,m 



l , m 
--In— + 2, 



5 



9m 2 

g 
(9m 2 



p (V ; o 

Pi (p 2 ; 0, m 



, m I \p2 = _ m 2 



1 /l , 
— — ~ + m — 

2m z \ e /i z 



giving the first term, 



Similarly we derive 



4 = ^ 2 



m z 



( 1 , m 2 \ 
- 3 (i + ln ^ +4 



/l , m 2 \ 
-3 ( | + ln-i ) i 



(40) 



(41) 



(42) 



Although not strictly necessary, we have carefully distinguished ultraviolet (1/e) poles from in- 
frared (1/e) ones, 

12 , 12 

— 7 — m 7T, - = — + 7 + in 7T, 



£ e £ £ 

with regulators that satisfy e + e = 0. The last term is the sum of three contributions, 

Bfj = i it 2 {8 e&j pi ■ pj C -A P_y • C u - C 12 P +ij - 2 B (P 2 ^.; m i? mj J, 

where the C-functions have arguments 

Pi = -eiPi, p 2 = P+ij, mi = 0,m 2 = mi,m 3 = m j; 

and where we have introduced 

P±ij = tiPi ± 
The function Co has a well-known representation, 

Co — „ fPi - + Fi 



where Pi j2 are given in terms of variables 

1 



ij 

2/1,2 



P 



P^ + m 2 - m 2 ± A 1/2 (-P 



2 2 



■5) 



Through the paper 



X(x,y,z)=x +y + z — 2(xy + xz + yz) 



(43) 
(44) 

(45) 
(46) 

(47) 

(48) 
(49) 
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represents the Kallen's A-function. The following result holds: 

1 



Fx 



p\n (v? - vf) 

Ph.- is 



In [ 1 — -4 I - In [ 1 - 4r 



F 2 = Fi In - ^ 
1 



+ 



1 



y\ J , 



P\ l3 (y\ j - yf) 



In 1 - 



In 



y, 



y(y - l)(x - yf 



f (yi,y l 2 J ) - f (y % 2,yi) 
'i-y 



-Li 



Lie 



y 



x - y 



x - y 



(50) 



Li 2 (z) is the standard di- logarithm. The higher rank three-point functions can be reduced to 
scalar ones as follows: 



- Cn €ipf + Ci2 P+jj = Cn + C Qpf + C12 ejpf. 
In the above equation the C-functions have arguments 

Pi = tiPi, P2 = CjPj, m 1 = TO;, m 2 = 0, m 3 = m r 
Following standard reduction techniques we obtain 

C U = - ^ ( m i Pi + e i € j Pi ■ Pj , 
1 

C12 = - - (m- P 2 + titj Pi ■ Pj Ri) , 



Pi = 



P. 



d 

Bo [Plij\ m h m j) ~ B o (- m f, 0, rrij) + /1 C ] , 

[-B (P+y-; mj) + P (-ml, 0, m,) + / 2 C ] , 
/1 = 2 m- , / 2 = - 2 e^- p* • p^-, d = ra-m 2 - (p 4 • p,)' 

giving the final result for P 3 , 



5J. = 2<^ 



4 e^- ^ • pj- C + P (P+y-; -J2 B o (-mf; 0, mj) . 



l=i 



If we use the explicit expressions for the scalar integrals we obtain 

1 



c = UU ln ^il 



2 \e 



Pi + 2*2 > 



(51) 



(52) 



(53) 



(54) 



P (-m 2 ; 0, m) = -I - In ^ + 2, P (P* - m,) = ~ - In ^ + P 3 , (55) 
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where P| est is the part of F 2 not proportional to Pi, see Eq. (|50|) . Furthermore, F 3 is the finite part 
of the .Bo-function, given by 



P 3 (p 2 ; 



rrii, rrij 



m 2 - m 2 m 2 A p 2 - i e + m 2 + m 2 - A 
— In — + 2 - — In J - , 

lp A rrij p z zmirrij 



with A 2 = A (— p 2 , m 2 , m 2 ^ . Using these results we find 

1 2 

B u = 2 iix 2 {2 (-1 + Pi ■ Pj - - J2 ln ^ 

l=*3 ^ 

+ F 3 }, 



+ 2 tit j pi - pj Fi ln 2 

A* 



P 2 ■ ■ - % £ 



(56) 



(57) 



showing that the virtual infrared pole originates from a Cq function. i*\ and F™^ are defined in 
Eq.flSTJ), P 3 in Eq.(|5§). Let Q^- = P+ij, we will consider two limiting cases. 



4.1 The case of large invariant 

If rrii = rrij = m the results simplify into 

2 



Ft 



nrcst 



In^ti 

P2.fl.. /?.. _ I' 



1 



In 



Pi 



2 £ Li 2 



'J 



i=±l 



2/% 



(5? 



with = 1 + 4m 2 /P 2 i7 -. Furthermore 



P 3 = 2 - A, ln 



Pij + 1 



In the limit | Qij 1^ m2 > where Qij = P+ij, we obtain 



Pi ~ T^r In 



j-ircst 
^2 



Pq ~ — 



m 2 



+ 



7T _ 



+ 2. 



(59) 



(60) 



13 



4.2 The case of small invariant 



Consider now the opposite limit where | Qfj |<^C m 2 ,m 2 . We easily derive 



Fx 



1 , m 2 

lri -L 

9 9 9 ' 

rrij — mi mi 



2 m 2 — mf 



Note that, for = m~ = m this further simplifies into 



1 9 ' 



1 



F 2 ~ — In—. 

m 2 /x^ 



Furthermore the two-point function is 
B (0;mi,mj) = ~ 
which, for equal masses, gives 



m 2 — mf 



2 .. m ) 9 . mf 



m„- In 



m s In ■ 



1 m 2 
- — In 



B (0;m,m) - ,, . 

e n z 

Collecting the various terms and introducing r = m 2 /m 2 we obtain 

-r + 1 



Bij = iir z 



r — 1 



lnr — 2 



1 r + 1 

e r — 1 



m 2 4-ln 2 4 



r — 1 /i 2 
which, for equal masses gives 



3 - r mf 1 - 3r , m 2 

ln^ + ln-f + 2f, 

r — 1 J 



B^ [Q 2 , mi = mj = m) — > 0, for Q 2 — > 0. 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



5 The two-particle radiation factor 

The photon spectral function is defined in Eq. (|2"2"D . Remarkably enough the exponent in E{x) 
can be computed exactly |pT[ . We will show it by computing E(x) for the case of only two 
emitters. This underlying ingredient will be referred to as the two-particle radiation factor, IZij. 
In dimensional regularization it reads as follows: 



1li 3 =^- n J d n k exp(-tk-x) 



5 + {k 2 ^ 



Pi ■ kpj ■ k ' 
where /i is the arbitrary unit of mass and 

P 2 i=-m 2 , q 2 j = (Pi + Pjf ■ 



(67) 



(6* 
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Furthermore 5 + selects positive energies, 

5 + (k 2 ) =6(k )5(k 2 ). (69) 
In order to evaluate TZij we introduce a parameter p^- defined by the following relations [Jl4f : 

p = pijPh q = Pj, (p-qf = Q- (70) 

Hence we have a solution 

1 



Pij 



2 m 2 



-q% - m 2 - m 2 ± \ 1 ' 2 (-gj, ro?> m 2 )} . (71) 



If q 2 j < we select p + , while p is chosen when q 2 j > 0. Next we introduce a Feynman parameter 
u and define 

P = q -\- {p — q) u. (72) 
With our choice for p it follows that P 2 < and we will compute 7£ in the frame where 

p = 0, P = M. (73) 

Therefore we obtain 



TZ — p duJ, 
Jo 



J= M 2 r(,/2-r) 7o 7-i M 1 "^) exp{-^r-* )}, 



M 2 T (n/2 - 1) Jo 



where r = Ixl. If we write 



then the function H becomes 

= y dkk n b J dy exp{- i k (yr - x )} [I - y ) 

~2\ n/2-3/2 



where J„ is a Bessel-function. As a consequence we arrive at the following result: 

M 2 ' 



ft = 2 "/ 2 - 3 / 2 p7T™/ 2 - 1 / 2 (pr) 4 "™ 



o 



(74) 



7r"/ 2_1 u 4_n r 1 dv 



i^-j ^ 2 V{n/2-l) J dkexp(ikxo) k^ 7 ' 2 J n/2 ^ /2 (kr), (76) 



K = K v + i e(x ) Ki = J™ dk exp (i ^fc) k n/2 ^ 7/2 J n/2 ^ 3/2 (k). (77) 
Real and imaginary parts of TZ will be computed separately. 
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5.1 The real part 

We start by computing the real part of K, Eq. ([77]) , 



dk COS^k) k n/2 ~ 7/2 J n /2-3/2(k), £ 



Xq 

r 



< £ < oo. 



(78 



This function will be considered separately in the two regions < £ < 1 and 1 < £ < oo. For the 
former we find 



o<e<i, K ^ = 2n/2 ~ 7/2 T r / (3/2) 2) 2Fl 



n-4 1 1 



2 ' 2'2' 



(79) 



where 2F1 is the standard hypergeometric function [|K]]. Let n = 4 + e', with e' > 0, then we can 
use the Laurent's expansion of the T-function 



£ 



--J + 0(s') 



50) 



and one of the transformation properties of the hypergeometric function to arrive at the following 
form: 



l-e'/2 



,2 2 2 

Since we are in the region < £ 2 < 1 the following relation holds: 



2 ' '2' ? , 



2^1 



2 ' '2' ? 



r(i/2) ~ r(f + i/2-eV2)r(z + i) g 
r (1/2-^/2) ^ r(z + i/2) 



/! 



(82) 



The function T (e'/2) shows an infrared pole and, therefore, we must expand the hypergeometric 
function in powers of e' . The whole procedure is cumbersome and, essentially, requires to obtain 
derivative of 2-^1 with respect to the parameters. Using the following expansion for the Euler 
gamma-function 

r (a + A e') = T (a) [l + A ^(a) e' + O (V 2 )] , (83) 
where it appears the Euler ^-function, we obtain 



E 

1=0 



r(2 + z) fl , e' 



/! 



{i + lka)-H^')]}^+°(^)- 



(84) 



As a consequence we are to consider the following expansion 



1F1 



l-e' 1 2 

"IT - ' 1; 2'* 



-1 £ 



(85) 



16 



The series in E 1 can be re-summed as follows. First we write 



and successively we obtain 



'dl!-'"— , 

1 — X 



l-x l 



^ 7o 1 - X 

1 — z Jo 



The re-summation gives a rather simple result, 



For K T in this region we find 



K r = (2tt)- 1/2 [- - 7 + ln2 - In (l - £ 2 ) - £ In - 



1 + C 



The result for 7?. r in the region < £ < 1 is 



= 71 P Jo I? " 7 + ln 77 + 2 ln 2 " ln ^ V ~ * ln I=f 



1 + C 



Now we turn to the complementary region 1 < £ < oo, where 



X r = 2 3 / 2 -«/ 2 £ 4 -" cos 



n - 4 \ r (n - 4) 



7T 



T(n- 1/2) 



7-2, 



71 — 3 71 — 1 



and where we will use 



2^1 



e' l + e' 3 + e' 



l-e'/2 



2' 2 ' 2 



l - 



2^1 



'3 3 + e' 
2' ' 2 



Since we are in the region where < £ 2 < 1 the following result holds: 



2^1 



3 3 + e 1 
2' ' 2 



l 



2 



m) + O (e' 2 ) , 



^(0 = g[^(i + |)-^(|)]e- 



2i 



The series for _E 2 can be re-summed by using the following identity: 



i)-1> 



1 , 1 — X 1 

dx x 1 ' 2 , 

o 1 — X 
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and, consequently, we get 



g [*('+§)-*(§ 



Jo 1 — X 



1=0 



dx 



,1/2 



1 — z Jo 1 — zx 



We obtain the following result for E 2 : 



and 



ft = [ I)" 2 [I - I 7 + I 1„2- i In (l - 1 ) - i{ lnl±i 



>7l7 Lg/ 2 2 

As far as the real part is concerned we have 

1 dit 



0<£< 1 ^ r = 7rp / — -- 7 + ln7r + 2 1n2-ln A tV-ein^-4 
jo M 2 l e' 1 — 4 



1 + e 



1 < £ < oo 1Z V = up 



r 1 du 
Jo M 2 



- - 7 + hivr + 2 ln2 - ln(-//V) - f In 



5.2 The imaginary part 

The imaginary part of the two-particle radiator is given in terms of 



Ki = J q °° dk sin k n/2 ~ 7/2 J n/2 ^/ 2 {k). 



This can again be written in terms of hypergeometric functions, 

0<£<1 K 1 = 2^I^y{^-) 2 P\(^=^,0;|;e 2 ), 

1< £ < oo ^ = 2 3 / 2 ""/ 2 , r( , W ~ 4) , , sin (n ~ 4)7F 

- s - s T (n/2 - 1/2) 2 



x 2 Fi 

We use the following expansion, 

(n - 4) 7r 



n — 3 n — 4 n — 1 



7T 



sin —— r (n - 4) = - T (n - 3) + O (n - 4) , 
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to show that, as expected, the imaginary part has no infrared poles and, therefore, we may set 
n = 4. The result is 



0<£<1 K l = 2-^^i 2 F 1 ^-\-e), 
It is straightforward to derive that 
The final result for the imaginary part reads as follows: 

°<e<i^ 2 /;^, 



(102) 



(103) 



(104) 



In order to cast the final result into a more compact form we introduce an infinitesimal quantity 
S such that 

x ^x + i5, 5 — > + , (105) 
It follows that x 2 — > r 2 — x\ — i x 5 and 

. , 2 \ i / 2 • r\ f ln(x 2 ) for x 2 > 

In(x -tin I — i^o = ^ i ) 2\ ■ I \r 2 r> 

v y V / [ ln(-a; 2 ) - «7re(a;o) for a; 2 < 

5.3 The complete result 

Collecting the results for 1Z = 1Z T + i e(x ) 1Z\, we obtain for < £ < 1 

Kr = pnj o ^[ i + 2(ln2- 7 )-ln / iV-£ln T -| 



^ 9 r 1 du _ 

= ^ Jo w ( - 



(106) 



and for 1 < £ < oo, 



1Z T = pit 



1 du 
M 2 

1 



Here we have introduced 



T + 2(ln2- 7 )-ln(-uV)-£ln . 
£ 4 - 1 J 



2 1 1 

— = - — 7 — U17T. 



(107) 
(108) 
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The total result is 

1Z = lZ r + i e(x ) IZi. 
For 1 < £ < oo it follows x 2 < and, therefore 

ln(— x 2 ) — i rre(x ) — > \nx 2 , 

while, for < £ < 1 we can replace 

In j— - + « 7re(x ) -> In ^— -y . 

Therefore the function 7?. is defined on the whole £-axis by the following expression: 

f 1 du 
71 = pit / — - 

Jo M 2 



(109) 
(110) 

(111) 



- + 2 (In 2 -7) -ln/iV -f ln^ii 
£ ? - 1 



P 
1 



•/ o 



' du r 2 2 £+ 1 

— [A IR -ln/xx -^ ln ^3T 



A IR = ^ + 2 (ln2- 7 ), 

e 

£ = — , x = a;o + i5, <5 — > 0_, 
r 



(112) 



To get the final form of our result we must express all quantities in covariant form. Therefore we 
have 



m 2 = -p 2 , x = -^f, e = 



(P-x) 



{P ■ x) z + P 2 x 2 



(113) 



6 Radiation from two legs 

Consider the case when the photon can be emitted by two external charged fermions only The 
corresponding radiator is 

9i6j 



2tt 2 



p 2 Tin + p 2 Kjj -2 pi- pj Kij 



and we need the photon spectral function, 

1 t 

Eij(x) = E(pi,pj]x) = -j——^ J d 4 K exp{-iK ■ x + alZipiiPj)}. 
We rewrite the radiator as 

n ij = Pij 71 ( A m - hi/iV) J + pijTrr-ij, 



= Plj 7t (A IR - \np 2 x 2 ) ^ In V — + Pij irr ij: 

1 p 2 

= - p tj 7T In (e- AlR p 2 x 2 ) -^—^ In — + Pij 7T r tj . 



(114) 



(115) 



(116) 
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In principle there is no problem in evaluating the finite part of the radiator explicitly, with a 
result that contains several di-logarithms. However, we have to exponentiate it as in Eq.( |115[ ) 
and, successively, we must compute the Fourier transform of the result: we will not be able to 
proceed any further with the complete expression. Using Eq. (|113|) we get 



and £ becomes 



P 



q 2 + (p 2 



q 2 ) u, 



P ■ x = q ■ x + (p — q) ■ x u, 



(117) 



(a + bu) 



Au 2 + 2Bu + C 

\2 



a = q ■ x, b = (p — q) ■ x, 



C = (q ■ x) z , B = q ■ x (p — q) ■ x + — (p 2 
Therefore, one integral is immediate 



2\ 2 
q ) x , 



.4 



1 du 
u M 2 



1 



p 2 — q 



1 P 2 

2 ln — » 



while the remaining one starts with 



du —h; ln — - - — 

u M 2 £ - 1 



du 



a + bu 



ln 



a + bu + U 



c + d-u) U a + bu — U 



118) 



(119) 



(120) 



where f/ 2 = lu 2 + 2Bu + C and, moreover, c = q 2 ,d = p 2 — q 2 . If needed the last integral can be 
computed through the substitution 

t = (121) 

U- 



C 



6.1 Coplanar approximation 

The Fourier transform of Eq.( |115| ) cannot be computed in closed form, therefore we change 
strategy and introduce an approximated formulas which is much simpler to handle in practical 
computations [|ll|]. This approximation is the coplanar one, where the effective photon momentum 
is constrained to lie in the plane formed by Pi and pj, so that the spectral function turns out to be 
proportional to S 2 (K^), where K± is the transverse component of K. This coplanar approximation 
reads 



^jln 

7T L 



_ g-A^ R A* Pi ' X Pj ' X 



1 + n.. K-Pi 

' ry 9 9 9 

Pj-PijPi 



P\ 



J 



where we have used the fact that for i 
introduced, 



3 Pi 



1,M 2 



2 6 

mf. Several new quantities have been 



A' 



IK 



1 3 
- - 27+ -, 

e ' 2' 



1 + 



Ql 



1/2 



miTrij , 



Q 



{eipi + ejpj) . 



(123) 
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Note that satisfies the following asymptotic behavior. 



Sij ~ 



Q 



1/2 







2 k s 2 2 
ij » //',.//',• 



The same result is rewritten as 

^ c (pup,) 
with a function Ajj defined by 



-Aij In 



.4; 



6i9j 

71 



pij 



Pi ■ Pj 



1 11 

mj - pfj mf 



In 



e -AJ? R A* Pi ' x Pj x 



Pij ml 



+ ^j, 



(124) 



(125) 




(126) 



There is some element of ambiguity in the definition of the coplanar factor of Eq. ( |122j ), which 
disappears when the difference (exact - coplanar) is properly included. It remains, however, when 
the result is expressed solely in terms of the coplanar approximation. This is connected to the 
fact that collinear logarithms, e.g. ln(Q 2 /m 2 ), are not fully accounted in the exponentiation and 
only double-logarithms, of the form \n(AE / E) ln(Q 2 /m 2 ), are properly included. 
In other words, what we can do is as follows: 



1. to exponentiate according to the YFS recipe and we do that, in principle, by including the 
full eikonal factor, no soft limit; 

2. to compensate with respect to the complete answer by including infrared safe residuals, (3% 
or more, see Eq. (|28|) ; 



3. to translate the bulk of exponentiation into structure functions (see Eq. ( |136|) below) including 
the rest into a remainder (see Eq.( |131| ) below) which could be handled numerically. 

From this point of view it really does not matter which approximation we start with, the only 
problem being that we, in general, do not control j3\ and neglect remainders. This is why the 
coplanar approximation is aimed to be as accurate as possible. 

The main characteristics that an approximation to the exact spectral function has to satisfy 
are: a) the possibility of extracting the typical form of the solution of the evolution equations 
for fermion(anti-fermion) distributions in the soft limit and b) the correct exponentiation of the 
leading logarithms. Furthermore, it should respect the correct scaling behavior and return no 
radiation for pi = pj. The particular choice made in Eq.( [L23| ) will be commented in Section || 

If we now write 



.4 



ij 



6i6j A. 



(127) 
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it is easily seen that Aij is non-positive Vi, j. Indeed we can immediately derive that 

•A-ij oc k /" dQfc *A.ij*A.^ij , 

^ = " ^e) • < 128 > 

where the integration is over the angular variables of the photon. From A^ • k — and A; 2 = it 
follows that .Ay • »4jj > 0. In coplanar approximation we have 

E c (p t , Pr ,K) = -^ J d A xexp{tK-x + an c ( Pl ,p,)}. (129) 

The total radiator will be the sum of its coplanar approximation and a remainder, 

n (pi, Pj ) = n c (pi, Pj ) + n icst (pi, Pj ) . (130) 

In this way we obtain the spectral function E(K) as 

Eij{K) = E (pi, P j] K) = j2~y J d 4 xex V {iK-x + an% PhPj )+aK™\ Pi , Pj )} 

= J d A K'<5>(K - K') j^-y J d 4 x exp{iK' -x + aTl }, 

/ d 4 yex P {iK-y + an rcst }. (131) 

The exact spectral function is now written as the convolution of some flux <3> with a kernel integral 
that we may cast in an appropriate form 



n(K) = -^-- J d 4 x exp{iK -x + a1l c ) 



(2 
(2 



l — A [ d 4 xexp(iK-x + a S tj ) \- f^PL^r£ ] ~ aA \ (132 ) 



The flux-function can be expanded in powers of a, giving 

*{K) = 5\K) + J d 4 y exp (iK-y) 7T cst + O (a 2 ) . (133) 

For the kernel Ti we use the fact that xq is defined with a small imaginary part, or 

ip ■ x — ■> ip ■ x + 5. (134) 

Futhermore, the relation 

I roo 

(ip ■ x)~ s = / da a 3 ' 1 exp (— ip ■ xa) , (135) 

T (s) Jo 
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is used to obtain 
H{K) = 



(2tt)* 

- A c 



d 4 x exp (i K ■ x + a 5i 



— e 



A?„ ^ Pi ' X Pj ' X 

1 H 



-a A; 



-aAi 



dado' (aa') aAij ~ 



T 2 (aAij 

Eg . (|136|) exhibits the typical form of a structure function, i.e. 

aA 



^3 

T (aA i:j + 1) 



„aAii — 1 



1 S A (a Pi + a' Pj - K) . (136) 



(137) 



Note that the quantity Aj R inside Eq.( |136|) exhibits the infrared pole. Eq.( |135| ) is valid for Res > 
and, therefore, Eg. ( |136j ) is valid only if KeAij > 0, which is not always the case due to the 6i6j 
factor. We will examine in the next section the generalization to positive exponents. 

The key relation, therefore, is Eq.( 135 ) which expresses the factor ip ■ x through the Mellin 
transform of an exponential that, in turn, re-establish energy-momentum conservation explicitly. 
For alternative uses of Mellin type representations see the first of ref . [|8] . 

To summarize our findings the coplanar approximation to the exact photon spectral function 
is given by 



E c { Pi , Pj ; K) = F m J d A K> j da t do 3 ${K' 
x 5 4 ioiPi + ajPj — K — K') , 

,2 



aAij 



r {aAij + 1 



\aAi 



e — 



-aA 



lJ e aSij . 



'138) 



There is an important property that any spectral function must satisfy, corresponding to the fact 
that a non-accelerated charged particle cannot radiate. From Eq . ([TT]) and from Eq. ( |128| ) it follows 
that both the exact and the coplanar spectral function satisfy 



(139) 



since both j M and A% are zero. 

A convenient inclusion of the flux-function $ towards the exact spectral function is achiev- 
able with the help of Eq.( |138| ): let us assume that both pi and pj denote outgoing momenta, 
then according to a well- know strategy which is adopted in the dipole formalism [15], we define 
transformed momenta pj,pj- Let us introduce the following quantities: 



A 



it 



\ ( ts2 2 2 2 2 
A ( — K , <7j 7«j , OaVTI 



A; 



A(- (a iPi + K' f ,-K 2 ,a 2 m 



(140) 



then the transformed momenta are 



1/2 



Pn 



Ay 

A,- 



Pn 



^Pl K 

K 2 P 



+ 2 



1 K 2 - a 2 m 2 + afm] 



K 2 



(141) 
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By construction they satisfy 

(TiPi + (TjPj + K' = pi + pj = K (142) 



and 



pi = -<r?m?, p) = -a)m). (143) 



Moreover, if the phase space element is written as 

d(j) (pi, Pj - K'\ K) = J] d* Pl 5 + (p* + mf) 6 4 (a iPl + a jPj - K - K') , (144) 

1=1,3 

after the transformation we get the following decomposition, 

d<P (pi,Pj] K'\ K) = d(j> {pi,pj- 0; K) ® [dK'), (145) 

which implies 

poo r- OlA ■ ■ n 2 

E (p^ Pj ; K) = F m / dadaj — - ,J (o-^-)^- 1 $ int 5 4 (Pi + ^ - ^) , 

JO L l (cuAy + 1) J 

$ int = | [^']$(^')- (146) 

The explicit form of the integration measure [dK'} depends on the nature of Pi,Pj, i.e. incoming 
and/or outgoing and the successive integration represents the difficult part of the procedure. We 
will not investigate it any further. 

Note that the coplanar approximation satisfies the scaling property 

E c ( Pi , Pj ; XK) = \^~*E C ( Pi , Pj ; K) , (147) 

and it is symmetric, i.e. E c (pi,pj] K) = E c (pj,pi] k). 

An additional comment concerns the correct interpretation of the parameters a, a' . In the 
standard approach the structure function represents the probability of finding an electron (or any 
fermion) within an electron (or any fermion) with a longitudinal momentum fraction a . In the 
standard approach, therefore, the emitted photons, typically representing initial state radiation 
in e + e~-annihilation are strictly collinear. There is a subtle point here [16]: consider e + e~ anni- 
hilation processes, dominated or not by s-channel diagrams. Let p be the four-momentum of the 
incoming electron in the laboratory system, 

1 77?^ 

p = -VS (0, 0, 0, 1) , /3 2 = l-4-^. (148) 

The electron, before interacting, emits soft and collinear photons. Let k — k\ + k-i + . . . be the 
total four-momentum of the radiated photons. Thus 

k=^{l-x) (0,0,1,1), (149) 
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so that k 2 = 0, as requested by collinear, massless, photons. Usually one can work with the massless 
approximation for the electron taking part in the hard scattering, thus an on-shell (massless) 
electron can emit a bunch of massless, collinear, photons and remain on its (massless) mass shell. 
But the electron mass cannot be neglected in particular cases and, after radiation, the electron 
finds itself in a virtual state having four-momentum 

p=p-k= -Vs (0,0,P-l + x,x), (150) 

with x being the fraction of energy remaining after radiation. As a consequence, the electron is 
put off its mass shell, 

p 2 = — ml + -(1 — j3)(l — x) s ~ — x ml for m e — > 0. (151) 

When considering the whole process we introduce p± for the incoming in the laboratory system. 
Once radiation has been emitted the momenta will be denoted by p± with 

P± = \V~s (0, 0, T(I3 - 1 + x ± ), x ± ) . (152) 
The total four-momentum becomes 

P = p+ + p- = l^y/s (0, 0, x„ - x + , x_ + x + ) , (153) 

with a corresponding invariant mass P 2 = —x + X- s = s. 

In our approach things are different since K inside Eq. ( 114(f ) is, by no means, restricted by 
the condition K 2 = 0. In other words K, instead of being collinear, is coplanar and we have no 
problem in dealing with relations like p = (1 — x)p and p 2 = —ml simultaneously. Needless to say 
coplanar includes collinear. 



7 Extension to n emitters 



The result of the previous sector should now be generalized to an arbitrary number of emitters. 
With n emitters we have n(n — l)/2 pairs and n(n — 1) a-variables. The corresponding spectral 
function becomes 

E(K) = — — J d A x exp (iK-x)H exp{aHf + aTZ\ est ), (154) 



i+i 



where = n(n — l)/2 and / runs over all pairings of external, charged, fermion lines. It follows 



E(K) = f d*K'§(K ! ){— A I d 4 x exp[i(#-JO + 



(2 7T) 



N 

d 4 y ex V (iK -y) Y[ex P {anr st }. 



(155) 



i=i 
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Consider now 

1 r N 

H N (K - K') = - - / (fx exp i (K - K') ■ x + a V Tl\ 

(2tt) 4 j l j=1 

1 ( \ ( 

= 77m / ^ Kn 777~74 / rf4x ex P( z ^ ' x + a ^} 

(Z7r) J (27TJ ^ 



/N— 1 
Aexp{i(ir-ir' - J fsT JV )- 2 / + a ^ ^}. (156) 



This result can be transformed into 

JV-l 



1 r t 

H N (K - K') = — - i / d 4 ^ / A exp{* (K — K' — K N ) ■ y + a ^ K?} 

roc ctA — 1 

x C N da N dcr' N ((r N (r' N ) a N 5 4 (a N p iN + cr' N p jN - K N ) . (157) 

•/ 

with 

c " = FW (158) 

This equation can be written as 

f°° aA 1 

Ti N (K - K') =C N dcr N da' N (a N a' N ) a N ~ H N -i (K — K' — a N p iN - a' N p jN ) , 

J 

H {K) =8\K). (159) 
The process can be iterated until we reach the final result 

H N (K -K') = Ha / daM M)^ 1 - 1 5 4 [K-K'-Y, + v'iPn) • (160) 
i=i Jo \ i=i I 

Define 

x k = ( a i + °i) k,jie{k}, 1 < k < n, (161) 
i 

where n is the number of emitters, N the corresponding number of pairs and {k} is the set of 
pairs that have one fc-line. Then 

N n 

Yl ( a iPk + a 'iPh) = J2 x k Pk- (162) 

1=1 k=l 

We recall that Eq. (|135|) is valid only for negative exponents, i.e. for positive Aij. Since some of 
the may be negative, depending on the product 9i0j, we are forced to consider an alternative 
derivation. Starting from 

E(K) = J d A K'${K')H{K - K'), 

H(K) — —^-4 f d 4 x exp{iK-x + aY, Ki}, (163) 
(2tt) J l 
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we derive 



H(K) = — U J] k A?R — 1 a ' / ^ exp (itf • x) ft (iPi ■ ^P Eie{i * } Al 
(2vr) l=1 % J i= i 



(164) 



As before /j denotes the set of all pairs that include the i-line. If the sum 

Ai= £ A t (165) 
ie{h} 

is positive the derivation follows as before. In the case of A; L negative we have to consider another 
integral representation, 



representation, 

^■ x) ~ aA ^V^A-)l^ da& 



.aAi-l 



exp (— iapi • x) — \ 



which is valid for 
By standard arguments it follows 

N 2 

H{K) = J] [e^ ^ 



1 < Re (aAi) < 0. 



(166) 
(167) 



x 1-9 I- £ A V ^ ) \ 5 * ^ - > 
where we have introduced a projector P 



Pfa) (j2 ~K^=5 4 (^2 a jPj - K^j 



(168) 



(169) 



7.1 The case 2^2 

For a 2 — > 2 process we have n = 4 external particles and iV = 6 pairs of emitters. Let us 
assume that all exponent are positive and consider 

(170) 



7ifi = 



6 roc ( N \ 

= UQ doM M^- 1 5*(k-k>-Y: (van, + a' lPjl ) 

1=1 J0 \ 1=1 ) 

With the following identification 



i — line j — line I — pair 

1 2 1 

1 3 2 

1 4 3 

2 3 4 

2 4 5 

3 4 6 



(171) 



28 



we obtain 



6 

(criPk + v'iPh) = (01 + o"2 + o"3) Pi + (el + 0-4 + 0-5) P2 

1=1 

+ (a' 2 +a' 4 + cr 6 ) p 3 + K + ^5 + O P4- ( 172 ) 



Therefore the object to compute is 



^ 6 = / n / n ^ 5 * - ^ - e 

i=i j ° 1=1 \ t=i 



XiPi 



X 8 (x± — (J\ — <7 2 - (T3) ^ (^2 — o{ — CT4 — Cr 5 ) 

x 5 (xg -a' 2 -a' A - (j 6 ) 5 (x 4 - 0-3 - 0-5 - o-g) . (173) 
With 1Z 6 expressed as the integral of some X 6 , 

^ = r n ^ 54 - ^ - e s e, (174) 
i=i \ i=i j 

we perform the various a integrations, starting with the trivial ones, 

r+00 6 6 

S 6 = / n da i II ( x i - °"2 - cr 3 ) aAl_1 (x 2 - cr 4 - cr 5 ) aAl_1 

J ° Z=2 i'=4 

x a^- 1 (x s ~a' A - af^ 1 (x 4 -a',- a',)^ 1 

x <t2 m ~ 1 (a' A ) aM ~ l °t Ah ~ X K)" 4 ^ 1 KT^ 1 • (175) 

After integration a set of consistency conditions will emerge, 

x 1 >a 2 + 0-3, x 2 > a A + <7 5 , 

^3 > ^4 + °"6, ^4 > cr' 5 + a' 6 . (176) 
The next step is to perform the ^-integration, 

f da 2 (a,)^- 1 (xl-a 3 -a 2 r A ^ 1 = B(aA 1 ,aA 2 ) ( Xl - o*)** 1 *^' 1 , (177) 

which requires x 1 > <7 3 and where B is the Euler's beta-function. Next we integrate over cr 3 , 

B(aA u aA 2 ) f° da 3 (a,)^' 1 (x, - a^^' 1 = ^ T ^ (^r^^ 1 ■ (i 78 ) 

Also the <t 4 , <7 5 integrations give 

/°° da 4 (a 4 ) oA4 - 1 (x 2 - <7 5 - a,)^ 1 " 1 = B (aAl, aA 4 ) (x 2 - a 5 ) a(Al+A4Kl , (179) 
Jo 
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B (aAi, aAi) I da^ (05) 
Jo 

r (aAi) T (al 4 ) r (aA 5 ) a{Al+Ai+A5 )-i 

\ X 2 



r (a {A x + 1 4 + A b )) 
Similar results hold for the <r 6 , a' 4 and a' 5 , o' % integration, with a total result 



180) 



n i= i r (aA {i} J i=1 
A {1} = A x + A 2 + A 3 , 1 {2} = A + 1 4 + As, 
A {3} = A 2 + A 4 + A 6 , 1 {4} =A 3 + A 5 + As, 



181) 



which can be easily generalized to 



N n ( \ , n 

E w = nr 2 (aAo [nr a £ aJ] _ n(^r Efce{fej 

i=l i=X \ ke{ki} J i=l 



182) 



where, for i fixed, the index k G {fcj} runs over all pairs that contain the line i and Vz,Zj > 0. 



8 Combining real and virtual corrections 

The emission of a real photon from the ij pair is described by a structure function language 
with an exponent a Ay — 1, with 



a _ Wj 

7T 



1 - Pij 



Pi ■ Pi ln P 2 iJ m h 



183) 



For virtual photons the overall, universal, factor 



B 



8tt 5 



K) ■ 



184) 



i<3 



is exponentiated and the relevant quantity is exp (2aKeB). For real radiation we have another 
overall exponentiation where the infrared-divergent object, for reach pair, is given in Eq.( |136| ), 



_ A c U 1-aA 

e m — 



1 e a5 ij 



exp | a Aij 



A^-lnf 



185) 



where / = {ij} is the emitting pair that we are considering. Here 



a; r = ,-2 7+ -, 



it 1 2 6 



186) 
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Structure function language means that the overall, real + virtual, exponent is multiplied by 

n mri- < 187 > 

where % runs over external charged lines. Moreover 

A = « E ( 188 ) 

where the sum is limited to those pairs containing %. 

8.1 IR finite exponent 

Once virtual and real exponentiation are combined we have cancellation of the infrared pole 
and some global remainder that reads as follows: 



a 



a 6, ; 6 



A^ + ln^" 


6i6j 1 







i a 



7T 



{ 



1 Pij o o 9 



J2 ^2 



In^i 



2 9i9j ReBij 

Sij \ 



A? R + ln^ 



(189) 



2 6 " J £ 

where 1/e shows the virtual infrared pole. The residue and the finite part of the virtual corrections 



are 



Bf= -l+aejPi-pjFl 3 , 



+ eiej pi ■ pj 



Ft In 



Q 2 -ie 
P 2 



+ K cst +- Fa- 



Here Fi is expressed as 



J/1,2 



2Q 



l - [Q 2 + mj - m? ± A 1 /' (-Q 2 ,m 2 ,m 2 ) 



F[ J = 



In (l - — ) - In ( 1- — 



Q 2 (yi - V2) 

Furthermore introduce a shorthand notation for the Kallen's function, 

A(-Q 2 ,m 2 ,m 2 )=A 2 . 



(190) 



(191) 



(192) 
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We see from Eq.(ffO|) that pij is also a solution of the equation 



If ej6j = +1 we find 



If instead e^j = — 1 we have 



Pij 



Pij 



2 mj 



(PijPi-Pjf = o. 



-Q 2 - m 2 - m 2 + A 



2 mi 



Q 2 + m 2 + m 2 + A 



In both cases we derive a noticeable relation, 

P 



mj — p 2 m 2 



Moreover, it is straightforward to show that 



1 



1 

T 
i 



Q 2 - 2/ 2 ) A' 



Consider now the quantity Y defined as 



Y 



2/i (2/2 



2/2 (2/i - 1) 
It follows that Y can be expressed as 

(Q 2 + m 2 + m 2 + A) 2 



Y 



4m 2 m 2 



4m 2 m 2 



(Q 2 + m 2 + m\ - AY 



Similarly we obtain 



m 2 



{Q 2 +m 2 i +m 2 j -Af . 



if €i€j = +1 



(Q2 +m 2 +rn 2 +A y _ 
-, 2 — T 11 c ? 'fc 9 ' — — 1 



In other words, an important result can be derived, namely 



ln^S = - e4 e,lnF, F - ^ - 1) 



m A 



2/2 (Vi ~ 1) ' 



The IR-finite exponent is therefore a/n 6i6j J 7 ^, with 

.2 1 



7T 



^• = l-2 7 -y + ln 



rriimj 



+ - ReF 3 + €i€ 



P i ■ Pj 
A 



In 



-- + 2 7 ] lnF + Re/ 2 rest }, 
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and with 

L = i n l&i F -t = 1 /r t (202) 



8.2 Asymptotic limits and general considerations 

a) the case Q 2 ^> m 2 

It is important to show the asymptotic behavior of this exponent in the region | Q 2 m 2 ,m 2 . 
We easily obtain that 

TYl 2 TYl 2 

A~-Q 2 , yi~-g|, (203) 
giving the asymptotic limit of F as 

2 2 

m-m- 

K (204) 

(g 2 ) 2 1 J 

Using the asymptotic behavior of Eq. (|60|) and also Eq. ( |124j ) we derive 



T l3 ~ 2 7 (L - - 1) - Z - Li j + 2 - ^ , (205) 

which shows, among other things, a cancellation of the In 2 terms. 

We can easily check that the exponent = a Aij — 1 has the usual asymptotic behavior. With 

s — — (pi + pj) 2 , rrii = rrij = m, Q$j = —1, (206) 

in the limit s 3> m 2 we get 

fln-^-lV (207) 
it \ m 2 J 

The result of Eq.( [205| ) explains the normalization in the definition of the coplanar factor. The 
infrared finite overall exponent has, in the asymptotic region | Q 2 m 2 , the correct behavior 
to reproduce the exponentiation commonly employed to describe initial state radiation in e + e~- 
annihilation, at least up to terms O (a 2 ) and without hard photons, 



G(x) = ^ x?- 1 exp{-/3 7 + 5 V+S }, 

a ( 3 , s 7r 2 \ 2 a f, s 



6 Y+S = - ir i n __ 2 + -l, (3 = — (In— - 1). (208) 
n \2 m 2 3 J it \ m 2 J 

In the above result G is the so-called radiator function which is connected to structure functions 
D by the relation 

G(x) = [ dzD(z)D(-). (209) 

Jx X 
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Note that Eq.( p08| ) is sometimes written as 

P 



G(x) = Yjf—J)^ 1 exp{-/? 7 + 5 YFS } (1 + 5 S ; 




*s = ^0+^ (-) In 2 —, (210) 
2 2 \7r J mi 



a form which follows from the evaluation of the YFS soft form-factor 

/?lne + 2 5 YFS , (211) 

e being the parameter introduced to limit the multiplicity of very soft photons. 

Moreover, if we neglect constant terms, our exponent of Eq. flZOlQ reproduces the leading be- 
havior of the Gribov-Lipatov solution ]T7|] of the evolution equation for the electron structure 
function, 

f^/2) 

which is valid in the soft limit. The factor Sij inserted in Eq. (|122|) has the purpose of reproducing 



D ^ = fmM\{l- 2 ^ "=- 61n ( 1 -^ to 4)- (212) 



the constant terms of Eq.( |208| ), therefore increasing the accuracy of the approximation. 

In conclusion the coplanar approximation, in the limit Xi -C 1, Vz coincides precisely with the 
exact expression resulting from the soft-photon re-summation, as given for e + e~ — > 7* in the 
classic YFS treatment. We want to stress that, for a general process, there are ambiguities in the 
choice of overall exponent in Eq. ( [20 1| ). Indeed the asymptotic factor 

+ (213) 

which is sub-leading, is tailored to reproduce the asymptotic form of the soft + virtual one loop 
corrections to the vertex e + e~7, i.e. 2 ReF dirac + <5 so ft- Single logarithms and constant terms in 
the overall exponent can never be exact, unless the full virtual + hard part of the spectrum is 
included. Therefore, the accuracy of the result is always limited to 

Px^ 1 [l + O (a L coll )}. (214) 

As mentioned in Section |2| the missing parts of the hard photon spectrum and of virtual corrections 
violate the well-known KLN result that the inclusive corrections are always small and free of large 
logarithms for a pair of final state emitters. Therefore, the accuracy of our result is, in this case, 
controlled only if tight cuts are imposed on the invariant masses of the final state pairs. To give a 
concrete example the cross-section for e + e~ — > ff that includes exact O (a) final state radiation 
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is fH 



a 



a c (s) = ^Qy(s)\-2(l- Z y + 4 



47T 



: — 2 In ( I - : ) I In — 



ni 



(215) 



+z 1 + - lnz + 2C(2) - 2Li 2 (1 - z) - 2 In (1 - z) + - - 3z 



where z = M 2 (ffj / s. The difference between the factor 



In (1 - z) 



In 



- 1 



m, 



(216) 



which is then exponentiated and the full result of Eq. (|216|) is due to hard, virtual and real, photons 
and is responsible for the correct limit z — > 0, 



Or 



a°(s) 1 



4 7T J 7 



(217) 



b) £/ie case Q 2 <C ?ri 2 
Finally we consider the case of one incoming/ one outgoing electron. For 



we obtain 



Pi 



1+ - 



t 



{Pi ~ Pit 

1/2 



-t <C mi 



mi 



a 



71 



1/2 



mi 



(2U 



(219) 



showing a power law behavior in the exponent, to be compared with the logarithmic one for 
-t > m\. 

When Q 2 <C m 2 ,m|, we derive 



A ~ ( 



2 2 \ 

—mA , 



Therefore the real emission is controlled by a coefficient 

,2 ,^2 



1 Pij 



Pi ■ Pj 



1 11 
171 3 ~ Pij m i 



ln Pii m i 



m 



1 + 2 



3 



m 2 + m 2 
l j_ 

1m 2 



(1 + r) 2 
4 - r 2 (l + r) 2 



In 



r(l + rl 



with r = m 2 /m 2 , which, for equal masses, reproduces the correct limit 



, mi = mj = m 



0, for Q 2 



0. 



(220) 



(221) 



(222) 



where a non-accelerated charge does not radiate. From these result we see that the normalization 
of the coplanar factor in Eq. (|122|) has been chosen to avoid appearance of spurious logarithms, 
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\n(Q 2 /m 2 ), for Q 2 <C m 2 . One may wonder whether the limit of very small momentum-transfer is 
relevant for any physical situation. Consider the process e + (p + ) + e~{pJ) — > e _ (g_) + X(q x ) and 
define 

Q = p--q-, y= E± ^- (223) 

P+-P- 

In the region of forward e~ scattering we have 



Q 2 > m 2 e ^— , (224) 



2 



where y is bounded by 

^<y<l-^, s = -( P+ +p_j\ (225) 



M 2 _ m e 



and Mo is the minimum invariant mass of the X-system. Whenever this mass is very low with 
respect to s we may reach the regime Q 2 <C m 2 .. 

As a final consideration, note that we are treating all charged fermionic lines on the same 
footing since only this combination of radiation is a meaningful gauge- invariant concept. For a 
general process, therefore, the interference between different legs is a fundamental part of the QED 
corrections and not an additional minor effect. Interference must be meaningfully definable, in 



particular when one exponentiates, see a discussion in ref. fL9fl . In certain situations interference 
is responsible for changing the scale in the leading logarithms but it should be present, as a 
matter of principle, also for those situations, as in forward scattering, where the typical scale is 
not large with respect to fermion masses. Strictly speaking the exponentiation of 'soft' initial-final 
interference is not accurate enough and, therefore, insufficient to describe annihilation processes 



around resonances. We will not dwell upon this subject any longer and refer to [20] 



9 Further refinements 

There are several reasons to increase the value of our approximation and ,for the sake of 
simplicity, we start our considerations by examining the case of a 2 — ► 2 process. The result may 
be cast into the following form, 

aocj dPS 2 J d 4 K$(K') J JJ dxt 

i=l 

x 5 4 (p+ +p--q+-q--K'-Y i x t p)j J] f(lT^j ^ (226) 
where (3o is the Born matrix element, = aAi and 

Pi = P+, P2 = P-, Ps = q-, Pa = q+. (227) 
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Furthermore (IPS2 is the two-body phase-space. If we neglect terms of O (a) in the flux-function, 
i.e. = 8 (K'), the argument of the delta-function in Eq. ( E2B ) becomes 

(1 - an) p+ + (1 - x 2 ) V- - (1 + a*) ?--(! + 2*) S- (228) 

Next we introduce scaled momenta, 

p+ = (l-x 1 )p+, v- = (1 - x 2 ) p~, 

<?+ = (1 + x 4 ) q+, q- = (l+ x 3 ) g_, (229) 

and derive four-momentum conservation in terms of the radiative process which incorporates the 
photons emitted along the directions of the charged fermions. There seems to be a clash between 
kinematics and matrix element; the original procedure requires a reorganization of the perturbative 
expansion which starts with the matrix element in soft approximation while the delta-function 
expressing conservation is kept exact, transformed with other ingredients into the photon spectral 
function which is again approximated to introduce conservation at the level of scaled momenta, as 
it would appear in the structure function language. In any intermediate step we are not authorized 
to use energy-momentum conservation since there is no delta-function to use. 

The lowest order factor in the perturbative expansion of the squared matrix element, (3q, is 
however constructed with non-scaled momenta. Typically, we will have 

f3 ccq^-p^q + -p + + q--p + q + -p_. (230) 

Let us change variables according to 

X + = 1 — Xi, X- = 1 — X2, 

— = 1 + 2:4, — = 1+0:3, (231) 

y+ y- 

so that the set of momenta satisfying conservation is 

p+=x + p +J p_ = x_p_, 

q+ - q- 
q+ = — , q- = —, 

y+ y- 

4 1 

TT dxi = 5 9 dx + dx_dy + dy_. (232) 

f=i y+y- 

As a consequence the Born matrix element can be cast into the following form: 

/?o oc [p_ ■q-p+-q++p+- q- P- • q+] • (233) 

x + x^ 1 J 

Futhermore,/?o will contain an overall factor s 3 from the s-channel propagator and from the root 
of the Kallen function. Given the relation between s and s, 

s = -( P+ +P-) 2 = —, (234) 
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we obtain 

Po = hI dipr ' (235) 

and, after transforming variables, 

J] / dx { p™ = / dx + dx-dy + dy_^— —^. (236) 
i=i s 

In deriving this result we have used scaled invariants, defined by 

s = ~ (p+ + P-f = -(q+ + q-f , 
i = ~ (p+ - q+f = ~ (P- ~ q-f , 

u = - (p+ - g-) 2 = - (q+ -P-f- (237) 

From phase space considerations and from positivity for all Xi it follows that 

< x±, y± < 1. (238) 

As a next step, we will show that it is possible to introduce, in a way that is consistent with the 
perturbative approach, a new lowest order result: it contains (3™ instead of (3™, 

f 2 _|_ f, 2 

= -p-, (239) 

the difference between the two formulations being of order a. This difference is known and 
computable so that perturbation theory indeed starts with a radiative kernel and the re-summation 
of all photons emitted along the directions of the external, charged, fermions. Let us consider this 
difference in more detail. First of all it is zero for x + = . . . y_ = 1. Next consider the following 
function: 

F(z,(3)= f 1 dxVp(x) f(x), (240) 

J z 

where (3 = a A — 1 and where we have also introduced the distribution 

V P = j3 (1 - x?- 1 . (241) 
Adding and subtracting a term we derive 

F (z, P) = /(l) (1 -zf + P I' dx f{x) - (1 - xf . (242) 

Jz 1 — X 

From this result the distribution can be computed. For instance, to second order in /3, we have 

V f) = 6(x-l) + (3{\n(l-z) + [-^—\} 

+ (? {\ In 2 (1 - z) + [^^] + } + O (f) . (243) 
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The '+ '-distribution is defined, as usual, by its action on a generic test function g{x): 

I dxg{x)f + {x) = f dx \g(x) - #(1)1 f(x). 

J z J z 1 

Consider now a simple example, where f(x) = x 2 f(x); it follows that 
F (z, f3) = f 1 (3 (1 - x f" 1 x 2 f{x) = f 1 dx(3 (1 - x f- 1 f(x) 

J z J z 

x {l + (3 [jL^ ^ dy (1 + y) f(y) + (l - x 2 ) In (1 - x)] + O ((3 2 ) }. 



(244) 



(245) 



Therefore, the perturbative expansion is controlled by the parameter f3 = a A — 1 and we may 
compute the kernel cross-section with scaled momenta and fold it with the appropriate factors of 
(3(1 — x) 13 ^ 1 . The difference with a kernel cross-section computed with non-scaled momenta and 
the successive application of the correct four-momentum conservation appears only at the next 
order in (3 and can be re-adjusted order-by-order in perturbation theory. 

There are two reasons why one should rescale momenta in the kernel cross-section. In any 
process with a resonance in the annihilation channel this procedure includes the possibility of a 
correct description of the radiative return directly in lowest order. There is more, this procedure 
is sometimes requested by a correct treatment of gauge invariance. Consider, for instance, the 
process e(p) + P(P) — > e(p') + X. We write 



da 



2P-p (,f) 



(2 7T 



,5+(p' 2 + m 2 e )5\p-q-p>) 



(246) 



The factor 5 4 (p — q — p') is successively promoted to become the spectral function E(p — q — p') 
which, in coplanar approximation, generates 



E c (p-q-p')^6 4 ((l-a)p-(l + a')p'-q), 
so that the kinematics is specified by 

P' 

q=p-p', p = xp = (1 - a)p, p' = - = (1 + a')p', 

y 



(247) 



(248) 



while the leptonic tensor in soft approximation is extracted as 



Tp, = \q 2 <V + p^p' nu + PuP 1 ^ = \q 2 V 



- T^u 4 ( — 



X 



2 



y 

X 



PuPnu + PvP'p, 



P^Pnu + PvP^ 



(249) 



and gauge invariance is respected only by T, namely q^T^ = q u T^ u = 0. 
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One should remember that, at least in principle, some of the exponents fti inside Eq. (|226D 
could be negative. In this case we should write 



n 



X 



Pi-l 



l-9{-^V{xi) fa 



(250) 



which, by setting some of the X{ to zero in the integrand, introduces the opportune subtraction 
on the original result. After changing variables, we obtain 



/ dx + dx_dy + dy_ (1 — x H 

Jo 1 ' 



i 



#4-1 



- 1 



X 



2 2 
X,X 



y+y- 



i-e{-(3 1 )v{ x+ -i) ... i-e(-fa)v{ y+ -i) A 



3inv 



(251) 



Each subtraction sets one parameter to one and corresponds, figuratively speaking, to disallow 
radiation from the corresponding leg. 

Furthermore, one should also take into proper account the kinematical constraints on the 
process. We must require that the invariant mass of the outgoing fermion - anti-fermion pair be 
at least 2 m f which implies 



x + x^y + y_ > 4 



m) 



(p++P-f 



(252) 



and x±,y± cannot be zero for non-zero fermion masses. 



10 A strategy for computing f3i 

The real improvement upon the present implementation of QED radiative corrections in generic 
2 — > n processes requires to go beyond /3 in Eq. (p8|) . Therefore, to go beyond the present 
approximation one has to compute f3\ in Eq. (|28|) or, at least to include the collinear singularity of 
the hard photon. Consider once more the process 

e + e- -vn/ + 7, (253) 

and let 

M = M^(k), M = M (e + e~ -> nf) . (254) 
Let % be an external fermion, for instance outgoing. Then 

M„ = teu( Pl ) (T; + 4) , T; = -i -Q- 7m fa + It) T(p t + k), (255) 
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where, for simplicity we have assumed massless fermions. T represents the contribution where 
the photon is emitted by the pj-hne with residual amplitude T l {j>i + k) and R represents the rest. 
Gauge invariance requires k ■ M = or 

ufa) k-tf = iQi ufa) T fa + k) . (256) 

Consider vectors Q, n and k±, with Q 2 = n 2 = and k± ■ Q = k± ■ n = and introduce 

k i Tin. 



kfi — z + k±^ 



2 z Q ■ n 



k 2 , n 



giving 

Using the relation 
we derive 



^Hl-)^-^-^^, (257) 

p 2 = k 2 = 0, 2 Pi -k = -—^-. (258) 
p i + k = Q + (k 2 ± ) , (259) 



£ \ M-e\ 2 =2Q 2 e 2 (1-z) 1 + z) * g | M ( Pi - + fe) | 2 +0 (1) . (260) 

spins -L spins 

The above result shows factorization of the collinear divergence. This representation solves many 
problems, for the leading k± behavior we must consider only external fermions and we do not need 
to have a precise knowledge of the residuals R % . Therefore we do not care about including internal 
lU-bosons emitting photons. 

However, the procedure is not gauge invariant, gauge violation occurring at O (/cjj, sub-leading 
w.r.t. leading In k 1 - corrections. There are two possibilities. Either photons are allowed only within 
a cone (with half-opening 5) surrounding each charged external fermion and we identify a leading, 
gauge-invariant, \n(E8/m) behavior with sub-leading gauge non-invariant contributions heavily 
suppressed or we integrate over the whole phase space of the photon. For the latter we may 
still identify and compute collinear logarithms without having to compute the exact O (a) matrix 
element but the scale in the logarithm becomes arbitrary. 

Therefore, a rigorous result cannot do without the exact, O (a) matrix element and any ap- 
proximation is not free from ambiguities. 



11 Conclusions 

One of the main ingredients in all calculations aimed to a very accurate control of high- 
energy-physics observables is represented by the re-summation of large QED corrections. This 
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procedure is usually performed by introducing structure functions. The scale that controls the 
large logarithms to be re-summed, as well as the i^-factor which one introduces to increase the 
accuracy of the calculation are based on some algorithm where one starts from the evolution 
equation for the structure function itself and seek for a solution which factorizes the re-summed 
Gribov-Lipatov [17] term and which, through an iterative method, matches with the finite, second 



order, result of Barbieri, Mignaco and Remiddi pl| . The procedure is to some extent based on 
the soft limit. 

The are are two classes of problems when one wants to generalize this algorithm to more 
complex processes, with many fermions in the final state. First of all the choice of the scale 
describing the evolution of the structure function is ambiguous and it is not even clear that one 
can have a realistic description with just one scale in situations where the dominant contribution 
to the process is far from the asymptotic regime. Secondly, an exact fixed order calculation is 
generally missing for the process and not only the scale in the dominant logarithms is ambiguous 
but also one has no control over the sub-leading logarithms. 

From this point of view, all claims that are based on general arguments as factorization the- 
orems or renormalization group equation are usually void. The only safe and rigorous argument 
that one can apply is the exact soft-photon re-summation, as given in the classic YFS treatment. 

The YFS algorithm has been further developed and modelled for its use in MonteCarlo pro- 
gramming and we have no pretension to be adding any substantial improvement, as clearly stated 
in the introduction. At the same time we cannot offer any claim pointing to the complete solution 
of the problem of a very precise implementation of QED radiation for processes more complex 
than e + e~-annihilation into fermion-antifermion pairs. This simple statement should not be con- 
fused with a failure of the method. In Section |3.2| we have repeated the classical argument that 
perturbation theory and exponentiation can me made consistent with residuals that are infrared 
finite. The accuracy at stake is confined in those ingredients that are missing just because of some 
technical inadequacy in controlling the full content of higher orders. 

What we have done is an attempt to systematize all arguments and speculations that have 
appeared in recent times in the literature. We started with the well-known YFS re-summation 
procedure, expressed in the modern language of dimensional regularization for infrared divergences. 
Based on this result we have adopted a slightly different approach where we, nevertheless, pursue 
the simple picture in which the whole effect of soft-photon emission is described by a real-photon 
spectral weight function. However, in our approach we avoid the introduction of a cutoff. In this 
respect we follow an old proposal by Chahine [JXXJ] by introducing an approximation to the exact 
spectral function which retains the important properties and incorporates the expected peaking of 
the emitted photons along the direction of charged particles. Also in this case we have completely 
reformulated the algorithm in modern language. 

The result of our investigation allows to write a corrected cross-section where the kernel for 
the hard scattering is convoluted with generalized structure functions where each of them is no 
longer function of one scale. Each external, charged, fermion leg brings a factor x aA_1 where a 
is the fine-structure constant and A is a function which depends on the momenta of the charged 
particles. 
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A preliminar account of these results has been given in |22 
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13 Appendix A 

The evaluation of the two-particle radiation function is based on Fourier cosine(sine) transforms 



of Bessel functions 23 



g(y) — I dx cos(xy) x 2fl 1 J 2u (ax) 



o 



, , 2 2 ^ 1 or 2 " T (a + v) „ ( 1 y 2 ' 

0<y<a, g(y) = f ^ + - _ -- 2 F, [u + -; - 

, , (a/2) 2u y- 2 "- 2 ^ T (2v + 2u) cos(z/vr + mr) 
a < y < oo, g(y) = 



T(2z/+ 1) 

1 a T 

2 y 2 



x 2 F 1 [v + n,v + n + -;2v + l- — \, (261) 



which is valid for 

3 

- Rez/ < Re/2 < -. (262) 
g(y) = I dx sin(xy) x 2 ^ 1 J 2u (ax) , 
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x 2 F l 
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-i r U + ^ + /^) 1 p A , . 1 3 ^ 

y r(l + I/ _^V 2 ^2 ^2 ^ '2'a 2 , 



. av _ 2At r(2 I / + 2 M ) ^ 



a < y < oo, = ( - ) y v * - - ^ sin(i/vr + /ivr) 



y 2 



Another useful integral is [23 



J + l dxexp(izx) [}-x 2 ) U 1/2 = 2 V rr 1/2 T (v + ^ z^J^z), (265) 
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(263) 



which is valid for 

1 3 

- Reu - - < Re/i < -, a > 0. (264) 



Rev > — \- 



which is valid for 

Further, the following integral representation holds [5 

if) (s + a) - if) (s + /3) = / dxx s ~ x , 

Jo 1 — X 



Res > — Rea, —Re/3. 



We have used two different integral representations for the T-function 153], 



z- s T(s) 



dxx s e 



1 „— zx 



Rez > 0, Res > 0, 



and a second one due to Cauchy and Saalschutz [53 



z~ s T(s) 



doc x> 



m=0 



m! 



which is valid for 



Re^ > 0, — (n + 1) < Res < -n. 



(266) 



(267) 



(268) 



(269) 



(270) 
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